Electrokinetic-flow-induced viscous drag on a tethered DNA inside a nanopore 
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Recent work has shown that the resistive force arising from viscous effects within the pore region 
could explain [Ghosal, S. Phys. Rev. E. 71, 051904 (2006) & Phys. Rev. Lett. 98, 238104 (2007)] 
observed translocation times in certain experiments involving voltage driven translocations of DNA 
through nanopores. The electrokinetic flow inside the pore and the accompanying viscous effects 
also play a crucial role in the interpretation of experiments where the DNA is immobilized inside 
a nanopore [U. Keyser, ei al. Nature Physics 2, 473 (2006)]. In this paper the viscous force is 
explicitly calculated for a nanopore of cylindrical geometry. It is found that the reduction of the 
tether force due to viscous drag and due to charge reduction by Manning condensation are of similar 
size. The result is of importance in the interpretation of experimental data on tethered DNA. 

PACS numbers: 87.15.Tt, 87.14.Gg 



The interaction of charged polymers, such as DNA 
with nanometer sized natural and artificial pores have 
received considerable attention recently [11-01 ■ Such stud- 
ies are partly motivated by the desire to understand how 
polymers cross internal membranes of cells [8j . The pos- 
sibility of developing devices capable of detecting prop- 
erties of biopolymers at the single molecule level for ap- 
plications such as rapid DNA sequencing Q is also a 
motivating factor for such studies. 

In a recent paper Keyser et al. reported [lOf experi- 
mental measurements in which a single strand of dsDNA 
was immobilized while threaded inside a nanopore by the 
application of a pulling force to counteract the electrical 
force on the DNA. This was achieved by attaching one 
end of the DNA to a Streptavidin coated polystyrene 
bead and holding the bead in a laser optical trap. The 
displacement of the bead from its equilibrium position 
could be detected and used to measure the pulling force 
on the DNA. The measured value was found to be about 
75 % of the maximum electric force on the DNA within 
the pore based on its bare charge irrespective of the elec- 
trolyte (KCl) concentration. This pulling force is how- 
ever determined by a complex interplay between elec- 
tric forces and hydrodynamics, as noted by Keyser et 
al. [lOj- The point of this calculation on an idealized 
physical model is to understand the relative importance 
of hydrodynamics and the reduction of effective charge 
on the DNA due to Manning condensation in determin- 
ing the observed pulling force. Since the DNA as well as 
the internal walls of the pore are charged, the pore region 
has a cylindrically symmetric distribution of oppositely 
charged counter-ions. In the presence of a strong electric 
field an electroosmotic fiow [11] is therefore generated 
in this region that fiows in a direction opposite to the 
direction in which the DNA would move if it were not 
immobilized (Figure [T|) . This fiow produces a hydrody- 
namic drag on the DNA partially balancing the applied 
electrical force. In this paper, a simplified geometry of 
the pore region is used to calculate explicitly the viscous 
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FIG. 1: Sketch illustrating the tethered nanopore experiment 
with a cylindrical pore 



drag. It is shown that the drag is a significant fraction of 
the total force acting on the DNA and needs to be taken 
into account for a proper interpretation of experimental 
data on DNA nanopore interactions. 

A simplified model is adopted in which the nanopore 
is regarded as a cylinder of radius i? (5.0 nm) and length 
L (60 nm). The part of the DNA inside the nanopore 
is regarded as a uniformly charged cylinder of radius a 
(1.1 nm) along the axis of the pore. The DNA has a 
linear charge density A (2 electron charges every 0.34 nm 
- the distance between adjacent bases) and a lower "ef- 
fective" c harg e density of Ag = A/gs due to the Oosawa- 
Manning [121 Il3 | phenomenon of counter- ion condensa- 
tion on its surface. The factor qB is the Oosawa-Manning 
factor, it has the value of qs — 4.2 for an ideal model of 
an infinite line charge in an unbounded electrolyte. Re- 
ferring to the system sketched in Figure [T] the fiuid ve- 
locity in the pore is axially directed and is described by 



some function u{r) where r is the distance from the cen- 
tral axis. The electric potential is —EQZ + (j){r) where the 
first term is due to the externally applied axial electric 
field Eq along the pore (the z-direction). The functions 
u and (f> are governed by the Stokes equation for viscous 
flow (with zero pressure gradient and an electric body 
force term) and the Poisson equation of electrostatics re- 
spectively: 
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where e is the permittivity of the electrolyte and pe is the 
electric charge density due to ions. The classical bound- 
ary conditions of 'no slip' are assumed for the velocity: 



u{a) = u{R) = 0. 



(3) 



Eliminating pe from the pair of equations ([T]) and ([2]) and 
using ([3]) determines u(r) in terms of the potential 0(r): 



u(rj = 



(r) - 0(i?) + A0 



ln(r/i?) 



ln(a/i?) 



(4) 



where A0 = <p{R) — <P{a). The viscous force (along the 
z-axis) on the DNA is then 



Fy = 2TTaLpu'{a) = 2TTaeEoL 
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If S is the surface charge density on the channel wall, 
then by Gauss' law, 



2'Kae(j)'{a) = Ae 
t(t)'{R) = S 



(6) 
(7) 



Using the first of these equations to eliminate (j)' {a) and 
noting that the electrical force on the DNA, F^ ~ XcLEq, 
equation ([5]) may be written as 
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where Ft = — F„ — F^ is the tether force. 

In order to calculate the ratio Fy/Ff. from equation ([8|), 
all that remains to be done is to calculate the quantity 
A(/). If the potential 4> anywhere in the gap does not 
greatly exceed fc^T/e (about 30 mV at room tempera- 
ture; fcs is the Boltzmann factor, T the absolute tem- 
perature and e the magnitude of electric charge) then 
may be computed from the Debye-Hiickel model (Ad is 
the Debye length): 
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with the boundary conditions ([6]) and ([7]). The solution 
to that problem is 
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FIG. 2: Ratio of viscous to electric force on DNA {~Fy/Fe) 
plotted as a function of KCl concentration assuming a charge 
of 2 electrons per base pair on the DNA (Manning factor 
of unity) and a constant surface charge on the pore wall 
parametrized by the surface charge density S (in mC/m^). 



where the constants A and B may be compactly ex- 
pressed in terms of the following dimensionless variables 
a■^, = a/Xo, R* = R/Xd and S'* = {2iTaS)/Xe- Thus, 



A 



B = 



S,Ki{a,) + Ki{R, 



a,[IiiR,)Ki{a,) ^ h{a,)KiiR,)] 

SJija,) + IijR,) 
a4h{R,)Ki{a,) - h{a,)Ki{R,)] 



(11) 
(12) 



where Kn, In (n is a non- negative integer) are modified 
Bessel functions of integral order. The potential drop A(j) 
in equation (|H]) is now easily found from equation (jlOp 
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where A/q = /o(^*) - 4 (a*) and AKq == Ko{R.,) - 
^o(a*)- 

The result of the calculation as described above is 
shown in Figures [H The figure shows the magnitude 
of the viscous to the electric force {—Fy/Fg,) evaluated 
from equation (jS]) using the value of Ad appropriate for 
a symmetric binary electrolyte [111 ]. The effective charge 
of the DNA is assumed to be the same as the bare charge 
(Ae = A) of -2e per base pair and the applied Voltage is 
AV ~ —120 mV. The surface charge concentration on 
the substrate, S has been assumed independent of the 
KCl concentration. From measurements of conductance 
and streaming potentials it has been shown [14, 15] that 
in Si/Si02 nanopores S ^ —60 mC/m^ for KCl con- 
centrations greater than about 0.1 M. However, for low 
concentrations, the surface charge density drops substan- 
tially and needs to be calculated from a more elaborate 
model that takes into account the equilibrium of surface 



reactions at the interface. For very low concentrations, 
S w —4 mC/m^. Due to the approximate nature of our 
model it does not seem worthwhile to attempt to incor- 
porate the proper dependence of S on KCl concentration. 
Instead, it suffices to show how —Fy/Fe varies with salt 
concentration for several fixed values of S between —4 to 
—60 mC/m^ as shown in Figure [51 The important fea- 
ture that these curves illustrate is that —F^/Fe is essen- 
tially constant for most of the KCl concentration range 
at a value of around 0.7 - 0.8. Thus, the viscous force 
is not small, and furthermore, if the viscous force were 
neglected and the reduction in the electrical pulling force 
were attributed to counter-ion condensation, it would ap- 
pear that the DNA effective charge is lowered by about 
75 percent (just the right amount to lead one to conclude 
that charge reduction by the Manning factor of gs =4.2 
is being validated!). Figure [2] shows an increase in the 
hydrodynamic drag with increase in the magnitude of 5, 
because the surface charge on the pore walls enhance the 
electroosmotic flow due to the DNA charge. A model that 
properly accounts for the variation of S with KCl concen- 
tration is expected to follow the S — —A mC/m^ curve 
closely for low concentrations (below about 0.1 M) and 
asymptote to the S = —60 mC/m^ curve at high concen- 
trations. The distribution of counter ions in the calcula- 
tions presented here was treated by means of the equilib- 
rium Debye-Hiickel theory and one may question whether 
that corresponds to the experimental conditions. Taking 
the ratio -Fy/F^ « 0.75, equation ^ gives A0 « 80 
mV if for Ag one assumes the DNA bare charge of two 
electronic charges per base pair. If this is reduced by the 
Manning factor of qs = 4.2, then Ac/i f* 19 mV. Though 
the formal requirement for the Debye-Hiickel lineariza- 
tion is \(j)\ « kuT je ~ 33 mV, in practice the double 
layer structure calculated from the Debye-Hiickel theory 
does not deviate substantially from the more accurate 
Poisson-Boltzmann calculation as long as the maximum 
value of \(j)\ is not substantially larger than 2kBT/e « 66 
mV [Ig. Therefore the Debye-Hiickel theory certainly 
suffices for our present purpose. For the cylindrical ge- 
ometry considered here, the applied potential does not 
disturb the equilibrium Debye layer structure, since the 
applied field is always along the iso-concentration sur- 
faces of the ions. However, for the real nanopore, the 
applied electric field may have a radial component, and 



one may ask if this is strong enough to distort the equi- 
librium Debye layer. Since \AV\ — 120 mV and L = 60 
nm, this imposed field is E^xt — \AV\/L ~ 2 x 10® V/m. 
The radial field within the Debye layer may be estimated 
as Eint = |A0|/(i? - a) w 5 X 10® V/m with Manning 
condensation and Eint ~ 20 x 10® V/m with the DNA 
bare charge. Thus, though the distortion of the dou- 
ble layer can be neglected for the purpose of obtaining a 
rough estimate, it should be accounted for if one desires 
an accurate calculation of the viscous force. In order to 
do so, the cross-sectional shape of the nanopore must be 
known. 

Analysis of this simplified model suggests that a more 
careful modeling is needed in order to properly interpret 
the Keyser et al. [10| experiments. Such a model should 
account for hydrodynamic drag while taking into account 
the proper pore shape, the variation of substrate charge 
with KCl, possible departures of the equilibrium poten- 
tial from the Debye-Hiickel model and other relevant con- 
ditions of the experimental set up. Numerical simulation 
on a more elaborate model incorporating these details 
used in conjunction with the experimental data could 
provide a more complete picture of the effective charge 
on DNA inside a nanopore. One may be tempted to 
question the use of classical continuum hydrodynamics 
to flows on the nanometer scale. However, the classi- 
cal approach has already been shown to give results in 
reasonable agreement with experiments on DNA translo- 
cation through solid state nanopores [3, [l3l . Molecular 
dynamic simulations, such as those presented by Aksi- 
mentiev et al. [l^ could be used to further reflne these 
calculations and to show that the effects described here 
persist even if continuum hydrodynamics is replaced by a 
discrete molecular model. In situations where the length 
of the DNA polymer is much greater than the length of 
the nanopore, entropic forces due to random coiling of the 
polymer become significant. Such entropic effects have 
been considered by Muthukumar |19l - l2l| . In the limit 
R ^ a, the hydrodynamic friction with the pore walls 
becomes unimportant and the problem becomes one of 
determining the electric field that would immobilize a 
polyelectrolyte acted upon by given non-electrical forces 
in the presence Brownian fluctuations. This problem has 
been studied in its general form by Long et al. [22, 123[ . 
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